Tunneling conductance of graphene NIS junctions 
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We show that in contrast to conventional normal metal-insulator-superconductor (NIS) junctions, 
the tunneling conductance of a NIS junction in graphene is an oscillatory function of the effective 
barrier strength of the insulating region, in the limit of a thin barrier. The amplitude of these 
oscillations are maximum for aligned Fermi surfaces of the normal and superconducting regions and 
vanishes for large Fermi surface mismatch. The zero-bias tunneling conductance, in sharp contrast 
to its counterpart in conventional NIS junctions, becomes maximum for a finite barrier strength. 
We also suggest experiments to test these predictions. 
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Graphene, a two-dimensional single layer of graphite, 
has been recently fabricated by Novoselov et. al. [1(1 . This 
has provided an unique opportunity for experimental ob- 
servation of electronic properties of graphene which has 
attracted theoretical attention for several decades [2| . In 
graphene, the energy bands touch the Fermi energy at 
six discrete points at the edges of the hexagonal Bril- 
louin zone. Out of these six Fermi points, only two are 
inequivalent; they are commonly referred to as K and K' 
points [H. The quasiparticle excitations about these K 
and K' points obey linear Dirac-like energy dispersion. 
The presence of such Dirac-like quasiparticles is expected 
to lead to a number of unusual electronic properties in 
graphene including relativistic quantum hall effect with 
unusual structure of Hall plateaus Recently, exper- 
imental observation of the unusual plateau structure of 
the Hall conductivity has confirmed this theoretical pre- 
diction Further, as suggested in Ref. @],the presence 
of such quasiparticles in graphene provides us with an 
experimental test bed for Klein paradox Q 

Another, less obvious but nevertheless interesting, con- 
sequence of the existence Dirac-like quasiparticles can be 
understood by studying tunneling conductance of a nor- 
mal metal-superconductor (NS) interface of graphene Q. 
Graphene is not a natural superconductor. However, su- 
perconductivity can be induced in a graphene layer in 
the presence of a superconducting electrode near it via 
proximity effect |{J E3] ■ It has been recently predicted 
[8| that a graphene NS junction, due to the Dirac-like en- 
ergy spectrum of its quasiparticles, can exhibit specular 
Andreev reflection in contrast to the usual retro reflec- 
tion observed in conventional NS junctions 11,12]. Such 



specular Andreev reflection process leads to qualitatively 
different tunneling conductance curves compared to con- 
ventional NS junctions |8|. However, the effect of the 
presence of a potential barrier between the normal and 
superconducting regions in graphene NS junction has not 
been studied so far. 

In this letter, we study the tunneling conductance of 
a normal metal-insulator-superconductor (NIS) junction 



of graphene in the limit of thin barrier. We show that in 
contrast to the conventional NIS junctions, the tunneling 
conductance of a graphene NIS junction is an oscillatory 
function of the effective barrier strength. The ampli- 
tude of these oscillations is maximum for aligned Fermi 
surfaces of the normal and superconducting regions and 
vanishes for large Fermi surface mismatch. In particular, 
we point out that the zero-bias conductance, in complete 
contrast to its behavior in conventional NIS junctions, 
reaches its maximum value for a finite barrier strength. 
By using the fact that the effective barrier height can be 
tuned experimentally by changing a gate voltage [!, 0] , 
we suggest an experimental setup where these effects can 
be observed. We also point out that our analysis re- 
produces the results of previous work on graphene NS 
junctions as a special case @. 

We consider a NIS junction in a graphene sheet oc- 
cupying the xy plane with the normal region extending 
x = — oo to x = —d for all y. The region I, modeled by 
a barrier potential Vq, extends from x = —d to x = 
while the superconducting region occupies x > 0. Such 
a local barrier can be implemented by either using the 
electric field effect or local chemical doping H 0] . The 
region x > is to be kept close to an superconducting 
electrode so that superconductivity is induced in this re- 
gion via proximity effect [1, Q . In the rest of this work, 
we shall assume that the barrier region has sharp edges 
on both sides. This condition requires that d <C 2n/kF, 
where fop is the Fermi wave- vector for graphene, and can 
be realistically created in experiments [6(. The NIS junc- 
tion can then be described by the Dirac-Bogoliubov-de 
Gennes (DBdG) equations @| 
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Here , ^„ = (^ fl , fea, ^ fi , -^a) are the 4 component 
wavefunctions for the electron and hole spinors, the in- 
dex a denote K or K' for electron/holes near K and K' 
points, a takes values K'(K) for a = K(K'), Ep denote 
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the Fermi energy, A and B denote the two inequivalent 
sites in the hexagonal lattice of graphene, and the Hamil- 
tonian H a is given by 



-ihv F {(T x d x + sgn(a)a y dy) . 



(2) 



In Eq. [21 vf denotes the Fermi velocity of the quasiparti- 
cles in graphene and sgn(et) takes values ± for a = K(K'). 

The pair-potential A(r) in Eq. [JJ connects the electron 
and the hole spinors of opposite Dirac points. We have 
modeled the pair-potential as A(r) = A exp(i(j))9(x), 
where A and <fi are the amplitude and the phase of the 
induced superconducting order parameter respectively 
and 9 is the Heaviside step function. The potential U (r) 
gives the relative shift of Fermi energies in normal, insu- 
lating and superconducting regions of graphene and can 
be modeled as U(r) = -U 9(x)+V 9(-x)9(x+d). At this 
stage, we introduce the dimensionless barrier strength 



X = V d/Hv F , 



(3) 



which is going to play a key role in the subsequent anal- 
ysis. In particular, we define a thin barrier as one with 
Vq — » oo and d — > such that x remains finite. For the 
NS junction studied in Ref. 0, \ = 0. The gate potential 
Uq can be used to tune the Fermi surface mismatch be- 
tween the normal and the superconducting regions. No- 
tice that the mean-field conditions for superconductivity 
is satisfied as long as Ao <C (Uq +Ef); thus, in principle, 
for large Uq one can have regimes where Ao > Ef 0. 

Eq. [1] can be solved in a straightforward manner to 
yield the wavefunction tp in the normal, insulating and 
the superconducting regions. In the normal region, for 
electron and holes traveling the ±x direction with a 
transverse momentum k y = q and energy e, the wave- 
functions are given by 

= (l,±e ±M ,0,0)exp[ l (±fc„x + gy)], 
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sin (a) = 



(0,0,l,Te ±W )exp^(±fc;.T + gy)], 



hvpq 



sin(a') = 



hvpq 



(4) 



e + Ep 1 K ' e — Ef ' 

where for the electron wavefunctions k n = 
(e + Ef) cos(a)/hvF and a is the angle of incidence 
of the electron. Similarly for the hole wavefunctions, 
k' n = {e — Ef) cos{a')/hvF with angle of incidence a'. 
Note that for an Andreev process to take place, the 
maximum angle of incidence for an electron is given by 
a c = arcsin[|e - E F \ / (e + E F )] [I]. 

In the barrier region, one can similarly ob- 
tain ip^ = (l,±e ±i0 ,0,0) exp[i (±k b x + qy)} and 

ijjjf = (o,0,l,Te ±lS ') exp[i(±k' b x + qy)} for elec- 
tron and holes moving along ±x. Here the angle 
of incidence of the electron(hole) 9(9') is given by 
sin [9{9')\ = hv F q/[e+(-)(E F -V )] and k b (k' b ) - 
[e- (+){E F - V )] cos [9 (9' )} /hv F - Note that in the limit 
of thin barrier, 9, 9' — > and k b d, k' b d — > x- 



In the superconducting region, the BdG quasiparticles 
are mixtures of electron and holes. Consequently, the 
wavefunctions of the BdG quasiparticles moving along 
±x with transverse momenta q and energy e has the form 



Te ±i( 7 -/J) -i* Te i(±7- 



sin(7) 



x exp [i (±k s x + qy) 
Hv F q/(E F + U Q ), 



(5) 



where k s — \J [(Uo + Ef) /Hvf} 2 — q 2 and 7 is the an- 
gle of incidence for the quasiparticles. Here k _1 = 
(hvF) 2 k s / [(Uq + Ef) Ao sin(/3)] is the localization length 
and P is given by 

= cos-^e/Ao) if |e| < A , 

= -icoBh _1 (e/Ao) if |e| > A , (6) 

Note that for |e| > Ao, K becomes imaginary and the 
quasiparticles can propagate in the bulk of the supercon- 
ductor. 

Let us now consider an electron incident on the barrier 
from the normal side with an energy e and transverse 
momentum q. The wave functions in the normal, insu- 
lating and superconducting regions, taking into account 
both Andreev and normal reflection processes, can then 
be written as [13] 



,5' ' 



(7) 



where r and r a are the amplitudes of normal and An- 
dreev reflections respectively, t and t' are the amplitudes 
of electron-like and hole-like quasiparticles in the super- 
conducting region and p, q, m and n are the amplitudes 
of electron and holes in the barrier. These wavefunctions 
must satisfy the appropriate boundary conditions: 



N \x=-d 



(8) 



Notice that these boundary conditions, in contrast their 
counterparts in standard NIS interfaces, do not impose 
any constraint on derivative of the wavefunctions at the 
boundary. Thus the standard delta function potential 
approximation for thin barrier [l2| can not be taken the 
outset, but has to be taken at the end of the calculation. 

Using the boundary conditions (Eq. [8|), one can now 
solve for the coefficients r and ta in Eq. [71 After some 
straightforward but cumbersome algebra, we find that in 
the limit of thin barriers, the expressions for r t, t' , and 
ta depend on the dimensionless coefficient x as 
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FIG. 1: Plot of tunneling conductance of a NIS junction 
graphene as a function of bias voltage for different effective 
barrier strengths for Uo = and Ao = O.OIEf- Note that the 
curves for \ = (black line) and \ = 7r(pink circles) coincide 
reflecting 7r periodicity. 




FIG. 2: Plot of zero-bias tunneling conductance as a func- 
tion of effective barrier potential \ f° r Uo — and Ao = 
O.OIEf- The n periodic oscillatory behavior with maxima at 
X = (n+l/2)n is to be contrasted with the monotonous decay 
of G(eV = 0) with increasing \ m conventional NIS junctions 



where the parameters T and p can be expressed in terms 
of 7, (3, a, and a' (Eqs.[H[5j and [6]) as 

= e~' 7 - 77 = e~ la ' cos(x) - zsin(x) 
e %1 + rj ' cos(x) — ie~ %a ' sin(x) 

e -i(7-/3) _ p e i(7-/3) 

P = re-^ + etf • (10) 

The tunneling conductance of the NIS junction can now 
be expressed in terms of r and r a by [l2j 

G(eV) = Go r (l-\r\ 2 + \r A \ 2 ^^)co S (a)da, 
Jo V cos(a) J 

(11) 

where Go = ^e 2 N{eV)/h is the ballistic conductance 
of metallic graphene, eV denotes the bias voltage, and 
N(eV) = (Ef + e)w / (wHvp) denotes the number of avail- 
able channels for a graphene sample of width w. Note 
that for eV <C Ep, Gq is a constant. Eq. [TT]can be eval- 
uated numerically to yield the tunneling conductance of 
the NIS junction for arbitrary parameter values. 

Eqs. [9] and [10] represent the key result of this work. 
From these equations, we find that in contrast to conven- 
tional NIS junctions [13], both r and ta are oscillatory 
functions of the effective barrier potential x f 0T an U Q^- 
gle of incidence a < a c . Consequently, we expect G(eV) 
(Eq. [TTj) to demonstrate oscillatory behavior as a func- 
tion of x with a period ir. We also note that our work 
reproduces the results of Ref. @ as a special case when 
X = riir for any integer n. 

Let us now consider the regime where the Fermi sur- 
faces of the normal metal and the superconductor is 



aligned (Uq — 0) and Ao -C Ep. A plot of the tun- 
neling conductances as a function of the bias voltage eV 
for different barrier strength x, shown in Fig. [TJ confirms 
the 7r periodic oscillatory behavior. The oscillation am- 
plitude is maximum for zero-bias (eV = 0), as shown in 
Fig. [21 and vanishes at the gap edge (eV = Ao). From 
Fig. [TJ we find two noteworthy features. First, the tun- 
neling conductance at the gap edge reaches a value close 
to 2G*o independent of the barrier strength, and second, 
the subgap tunneling conductance becomes close to 2Gq 
when x = {n+ l/2)7r for any integer n and any eV < Ao. 

Both of the above-mentioned features can be under- 
stood by noting that when eV <C Ep, a ~ —a' ~ 7 
(Eq. 2]). In this limit, using Eqs. M and [TU1 it can be 
shown that the reflection amplitude, for eV < Ao, 
becomes r ~ M [cos (x) + i sin(x) cos(a)] /T>, where V = 
[cos(a) ]cos(/3) cos(2x) + i sm(f3)] + 2 sin(2x) sin(/3) sin 2 (a)] 
and 

2 sin(a) [sin (x + f3) — sin (x — (3)] . (12) 

Note that A/" and hence r vanishes and when a — or 
sin(x + P) = sin(x — (3). The former condition (a — 0) 
is a manifestation of the well-known Klein paradox [7| 
and occurs since scattering of normal-incident Dirac elec- 
trons from a potential barrier can not change their chi- 
rality This effect, however, is not manifested eas- 
ily in the tunneling conductance since G receives con- 
tribution from electron approaching the barrier with all 
possible incidence angles a < a c . The latter equality 
(sin(x + 0) = sin(x — 0) ) represents condition for trans- 
mission resonance (r = and \ta\ = 1) of a graphene NIS 
interface and has the solutions [3 = and x = (n+l/2)7r. 

Such transmission resonances occur for all barrier 
strengths x and angle of incidence a when eV — Aq {i.e. 
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FIG. 3: Tunneling conductance as a function of bias voltage 
for effective barrier strengths x = (solid line) and % — n/2 
(open circle) for Uo = 25Ef and Ao = 2Ef- The tunneling 
conductance becomes barrier independent and vanishes for 
eV — Ef- The inset shows a schematic experimental setup. 
The dashed region sees a variable gate (shown as pink filled 
region) voltage Vo which creates the barrier. Additional gate 
voltage Uo, which may be applied on the superconducting 
side, and the current source is not shown to avoid clutter. 

(3 = 0). Consequently, G(eV = Ao) ~ 2Gq independent 
of the barrier strength, as is also well-known for conven- 
tional NIS junctions [12]. The novel aspect of a graphene 
NIS junction comes from the second class of solution of 
the transmission resonance condition: \ = i n + V^) 71 "- 
At these special values of the barrier, G ~ 2Go for any 
subgap voltage as long as eV <C Ep [13| . This fea- 
ture, clearly seen in Figs. [T] and [21 is in sharp contrast 
to conventional NIS junction where G(eV < Ao) al- 
ways decreases with increasing barrier strength [13] ■ For 
X i= (n + l/2)7r or eV A , r ^ and consequently 
\ta\ < 1 so that G < 2Gq. In particular, r reaches a 
maxima leading to minimum value of subgap tunneling 
conductance for \ — n7r - Thus, in contrast to conven- 
tional NS junctions, G(eV = 0) < 2G for y = in 
graphene NS junctions as noted earlier in Ref. [8J. 

Next, we briefly explore the regime where Ao > Ep 
and U > E F , so that A < (U + E F ) [8]. Here the 
tunneling conductance, shown in Fig. [3] for Ao = 2Ep 
and x = 0, 7r/2, becomes independent of barrier strength. 
This can be intuitively understood from the fact that a 
large Fermi surface mismatch acts as an effective barrier 
for the electrons tunneling at the interface [2] which 
makes the presence of an additional barrier irrelevant. 
Further, as seen from Fig. O the tunneling conductance 
vanishes for eV = Ep = 0.5Ao due to the absence of 
Andreev reflection since a c = at this bias-voltage. Our 
results in this limit therefore becomes identical to those 
for NS junction studied in Ref. [§]. 



Finally, we discuss possible experimental setup (shown 
in inset of Fig. [3]) to test our predictions. In the exper- 
iment, the local barrier can be fabricated using meth- 
ods of Ref. [f|. The easiest experimentally achievable 
regime corresponds to Ao <C Ep with aligned Fermi sur- 
faces for the normal and superconducting regions. We 
suggest measurement of tunneling conductance curves at 
zero-bias (eV = 0) in this regime. Our prediction is 
that the zero-bias conductance will show an oscillatory 
behavior with change of effective bias voltage with max- 
ima(minima) when x becomes odd(even) integer multi- 
ples of 7r/2 . In graphene, typical Fermi energy is Ep ~ 
80meV and the Fermi- wavelength is A = 2ir/kp ~ lOOnm 

0] . For realization of the thin and sharp barriers dis- 
cussed in this work, one needs d/X « 1 and Vo/Ep 1. 
Effective barrier strengths of 500 — lOOOmeV and barrier 
widths of d ~ 20 — 10 nm, which can be achieved in real- 
istic experiments 0, @] , shall therefore meet the demands 
of the proposed experimental setup. To observe the oscil- 
latory behavior of the zero-bias tunneling conductance, 
it would be necessary to change \ in small steps 5\- For 
barriers of a fixed width, with values of d/X = 0.1 and 
Vo/Ep = 10, this would require changing Vo in steps of 
approximately 12meV which corresponds to 8\ = 0.1. 

In conclusion, we have a presented a theory of tunnel- 
ing conductance of graphene NIS junctions. We have 
demonstrated that the tunneling conductance exhibits 
novel 7r periodic oscillatory behavior as a function of bar- 
rier strength of the junction and have suggested exper- 
iments to observe this effect. The authors thank S.M. 
Bhattacharjee, A. Ghosh, P.K. Mohanty,M. Khan, T. 
Senthil and V. B. Shenoy for valuable discussions and 
Graduate Associateship Program at Saha Institute which 
made this work possible. 
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